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Abstract:  This  report  discusses  major  developments  that  have 
occurred  since  the  release  of  NSM  Versioni.o.  These  improvements 
were  in  the  areas  of:  l)  partitioning  of  nutrients  (dissolved  and 
adsorbed)  within  surface  runoff,  2)  mass  transfer  of  nutrients 
between  the  soils  and  overlying  water  column,  3)  erosion  and 
sedimentation  of  particulate  nutrients,  and  4)  soil  plant  dynamics. 

Nutrients  at  or  near  the  soil  surface  can  be  transformed  to  overland 
flow  in  solution  form  by  the  mixing  of  rainwater  with  soil  solution.  A 
rate-limited  mass  transfer  at  the  soil  overland  flow  interface  controls 
the  dissolved  nutrient  transfer  from  soil  solution  to  overland  flow, 
once  formed.  The  mass  transfer  coefficient  will  be  affected  by  many 
factors  including  characteristics  of  the  raindrops,  overland  flow,  and 
nutrient  characteristics.  The  predictions  of  these  variables,  together 
with  nutrient  transport  in  the  surface  runoff,  need  to  be  validated 
with  the  measured  data  at  the  field  scale. 

Regarding  soil  plant  dynamics  capabilities  within  NSM,  the  EPIC 
formulations  appear  to  provide  realistic  responses  for  simple  annual 
based  systems.  It  is  necessarily  simple  for  use  with  the  available 
datasets  and  should  be  improved  by  addition  of  plant  growth  pro¬ 
cesses  with  nutrient  cycling.  Further  development  of  this  module  for 
a  more  detailed  treatment  of  plant  water  uptake  and  associated 
nutrient  interactions  is  needed.  Attention  to  refinement  of  links 
between  soil  plant  dynamics  and  hydrological  processes  is  also 
required. 

As  research  continues,  improved  understanding  of  nutrient  process 
descriptions  and  new  capabilities  will  be  developed  and  integrated 
into  the  NSM.  Future  releases  of  NSM  will  include:  1)  improved 
nutrient  kinetics  for  channels,  2)  improved  soil  plant  dynamics 
(EDYS  sub-module),  3)  inclusion  of  carbon  cycle  modules  for  soils 
and  overland  flow,  and  4)  improved  nitrogen  and  phosphorus 
kinetics  for  soils  and  overland  flow. 
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Preface 

This  report  summarizes  efforts  undertaken  in  developing  a  nutrient  sub¬ 
model  for  linkage  with  a  variety  of  hydraulic  and  hydrologic  modeling 
systems  since  the  release  of  SWWRP-NSM  Version  l.o.  This  development 
effort  was  performed  by  the  Engineer  Research  and  Development  Center 
(ERDC),  Waterways  Experiment  Station,  Vicksburg,  Mississippi.  Funding 
was  provided  under  the  System-Wide  Water  Resources  Program 
(SWWRP).  Appreciation  is  extended  to  all  those  who  assisted  in  the 
formulation  and  review  of  the  process  descriptions  implemented  within 
the  Nutrient  Sub-Model  (NSM). 

Principal  Investigators  for  this  study  were  Dr.  Billy  E.  Johnson  of  the 
Water  Quality  and  Contaminant  Modeling  Branch  (WQCMB),  Environ¬ 
mental  Laboratory  (EL),  U.S.  Army  Engineer  Research  and  Development 
Center  (ERDC),  and  Terry  K.  Gerald  and  Dr.  Zhonglong  Zhang  of  SpecPro, 
Inc.  Gerald  and  Dr.  Zhang  were  funded  as  onsite  contractors  under  Task 
Order  Contract  W912HZ-05-D-0011  on  Civil  Delivery  Order  No.  15. 

Dr.  Johnson  conducted  his  portion  of  the  study  under  the  general  super¬ 
vision  of  Dr.  Barry  W.  Bunch,  Chief,  WQCMB;  and  under  the  general 
supervision  of  Dr.  Richard  E.  Price,  Chief,  Environmental  Processes  and 
Engineering  Division,  EL;  and  Dr.  Beth  C.  Fleming,  Director,  EL. 

COL  Gary  E.  Johnston  was  Commander  and  Executive  Director  of  ERDC. 
Dr.  James  R.  Houston  was  Director. 

This  report  should  be  cited  as  follows: 

Johnson,  B.  E.,  T.  K.  Gerald,  and  Z.  Zhang.  2008.  System-Wide 
Water  Resources  Program  Nutrient  Sub-Model  (SWWRP-NSM) 
Version  1.1.  ERDC/EL  TR-08-25.  Vicksburg,  MS:  U.S.  Army 
Engineer  Research  and  Development  Center. 
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1  Introduction 

Non-point  source  (NPS)  runoff  of  pollutants  is  viewed  as  one  of  the  most 
important  factors  causing  impaired  water  quality  in  freshwater  and 
estuarine  ecosystems  and  has  been  addressed  as  a  national  priority  since 
the  passage  of  the  Clean  Water  Act,  Section  319  (33USC  $1329).  Nutrients 
are  conveyed  to  water  bodies  from  watershed,  atmospheric,  and  oceanic 
sources.  To  comply  with  water  quality  regulations  as  well  as  long-term 
watershed  planning  and  management,  advanced  watershed  water  quality 
models  need  to  be  developed. 

Within  the  System-Wide  Water  Resources  Program  (SWWRP),  multiple 
riverine,  estuarine,  watershed,  and  subsurface  flow  models  are  being 
modified  to  address  issues  of  environmental  concern.  Several  integration 
approaches  are  either  ongoing  or  proposed  to  accomplish  this  task.  To 
have  a  full  system-wide  water  quality  and  contaminant  capability  in 
SWWRP,  the  different  hydrologic  and  hydraulic  engines  must  utilize  a 
common  water  quality  and  contaminant  approach  to  prevent  the  arbitrary 
partitioning  of  state  variables.  The  goal  of  this  development  effort  is  to 
upgrade  existing  hydrologic  and  hydraulic  models  (i.e.,  water  engines) 
using  a  common  water  quality  approach  in  order  to  facilitate  their  linkage 
and  application  on  a  system-wide  basis. 

In  keeping  with  a  common  water  quality  approach  to  model  development, 
a  library  of  water  quality  kinetic  modules  are  being  developed  such  that 
they  can  be  integrated  with  a  variety  of  water  transport  engines.  The 
library  of  algorithms  has  the  following  characteristics: 

•  Multi-species,  multi-phase,  and  multi-reaction  system 

•  Fast  (equilibrium-based)  and  slow  (non-equilibrium-based  or  rate- 
based)  reactions 

•  Easily  extensible  to  new  reaction  pathways 

•  Includes  both  common  nutrient  and  contaminant  packages  as  well  as 
geochemistry 

•  Simple,  well-defined  data  interface  and  calling  procedure. 
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The  water  quality  modules  are  developed  such  that  they  are  data  structure 
independent,  thus  facilitating  their  integration  into  a  wide  range  of 
modeling  systems. 

The  Nutrient  Sub-Model  (NSM)  has  been  developed  as  a  library  of  water 
quality  kinetic  modules.  NSM  considers  detailed  Nitrogen  (N)  and 
Phosphorus  (P)  cycling  and  computes  nutrient  kinetic  fluxes  for  nitrogen 
and  phosphorus  at  the  watershed  scale.  Modeling  of  nutrients  within  NSM 
consists  of  three  distinct  parts: 

•  Simulating  the  N  and  P  cycle  in  the  soil 

•  Transformation  and  loading  of  N  and  P  species  in  the  overland  flow 

•  Simulation  of  the  N  and  P  cycle  in  the  water  column  (both  overland 
and  channel). 

NSM  is  written  in  a  modular/process  structure  facilitating  the  coupling 
with  transport  components  from  a  hydrologic  and  hydrodynamic  (H&H) 
model  (e.g.  HEC-RAS,  ADH,  HEC-HMS,  GSSHA).  Concentrations  for  each 
nutrient  species  are  simulated  spatially,  and  are  updated  at  each  compu¬ 
tational  time  step,  along  with  water  flow  and  solute  transport. 

This  report  describes  partitioning  of  N  and  P  between  the  dissolved  state 
and  the  adsorbed  state  in  overland  and  stream  flow,  the  mass  transfer  of 
nutrients  between  the  soil  layer  (or  river/lake  bed  layer)  and  the  overlying 
water  column,  the  erosion  and  sedimentation  of  nutrients,  and  the  current 
soil  plant  dynamics  available  within  NSM.  Nutrient  kinetic  process  theory 
and  formulations  for  overland  and  channel  flow  are  described  in  detail  in 
Johnson  and  Gerald  (2006)  and  will  not  be  discussed  in  this  report. 
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2  Partitioning  of  Nutrients  in  Surface 
Runoff 

Nitrogen 

Soil  nitrogen  occurs  in  several  chemical  forms,  which  vary  greatly  not  only 
in  their  characteristics  but  also  in  behavior.  Soil  nitrogen  forms  are  classi¬ 
fied  as  either  organic  or  inorganic  nitrogen.  Typically,  most  nitrogen  in 
soils  and  surficial  sediments  occurs  in  organic  form.  Organic  forms  of 
nitrogen  are  found  in  compounds  such  as  amino  acids,  protein,  and  more 
resistant  nitrogen  compounds  (ultimately,  humus).  Organic  nitrogen  pools 
are  differentiated  on  the  basis  of  their  carbon-to-nitrogen  ratio  (C:N). 
Organic  nitrogen  and  NHA  are  mostly  absorbed  by  clays  and  can  be  moved 

with  the  soil  during  erosion.  In  such  forms,  nitrogen  can  be  considered 
immobile.  Inorganic  forms  include  ammonium  ( NHA  ),  nitrate  ( N03 ), 

and  lower  concentrations  of  nitrite  ( N02 ).  More  importantly,  however, 
NHa  and  N03  are  soluble  and  are  mobilized  through  the  soil  profile 
during  periods  of  rain  by  the  process  of  leaching. 

The  dominant  nitrogen  species  in  waters  are  dissolved  inorganic  nitrogen  - 
NHj ,  N02 ,  N03 ,  dissolved  organic  nitrogen,  particulate  organic  nitrogen 

and  particulate  inorganic  nitrogen  (Burt  et  al.  1993).  The  following  nitro¬ 
gen  state  variables  within  the  water  column  are  simulated  by  NSM: 

•  Particulate  Organic  Nitrogen 

•  Particulate  Ammonium  Nitrogen 

•  Ammonium  Nitrogen  ( NHA  ) 

•  Nitrate  Nitrogen  ( N03 ) 

Phosphorus 

Soil  phosphorus  can  exist  in  a  variety  of  forms:  organic  phosphorus  bound 
with  carbon  and  oxygen  in  plant  matter,  inorganic  phosphorus,  and  as 
dissolved  (soluble)  phosphorus.  The  major  types  of  organic  phosphorus 
compounds  are  added  to  the  soil  by  recycling  plant,  animal,  and  microbial 
remains  in  the  soil  environment.  Inorganic  phosphorus  occurs  as  primary 
minerals  (those  derived  directly  from  weathered  parent  material)  and 
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secondary  minerals  (those  formed  by  precipitation  of  phosphorus  with  Al, 
Ca,  and  Fe).  Organic  phosphorus  and  inorganic  phosphorus  are  associated 
with  the  particulate  phase.  In  studies  of  phosphorus  movement  from 
agricultural  lands  the  majority  is  sorbed  onto  clay  materials  and  trans¬ 
ported  as  erosion  products.  Soluble  phosphorus  exists  largely  as  ortho¬ 
phosphate  ions  and  can  vary  from  o.oi  to  3.0  mg/L  with  organic  phos¬ 
phorus  compounds  also  present  in  solution.  Phosphorus  is  generally  much 
less  mobile  than  nitrogen,  being  strongly  adsorbed  to  soil  and  organic 
matter.  Consequently,  in  aquatic  studies,  the  focus  is  often  on  the 
sediment-associated  forms  of  phosphorus,  as  these  tend  to  dominate  the 
total  phosphorus  flux. 

Phosphorus  enters  surface  water  primarily  as  particulate  matter  and 
secondarily  as  dissolved  inorganic  phosphorus,  also  known  as  “ortho-P” 
(H3PO4  and  its  conjugate  base  forms).  The  following  P  state  variables 
within  the  water  column  are  simulated  by  NSM: 

•  Particulate  Organic  Phosphorus 

•  Particulate  Inorganic  Phosphorus 

•  Dissolved  Inorganic  Phosphorus 

Partitioning  of  inorganic  nitrogen  and  phosphorus 

The  physical  processes  of  surface  runoff  cause  pollutants  or  natural  mate¬ 
rials  to  be  transported  and  mixed,  or  exchanged,  with  other  media.  There 
are  primarily  two  kinds  of  physical  processes  by  which  pollutants  are 
transported  in  fluids:  bulk  movement  of  fluids  from  one  location  to 
another  (advection),  and  random  mixing  processes  (diffusion)  within  the 
fluids.  This  random  diffusive  process  is  often  termed  Fickian  transport 
(Fischer  et  al.  1979).  Besides  dissolved  species,  pollutants  that  are 
absorbed  on  the  surfaces  of  particles  or  absorbed  into  particles  can  also  be 
affected  by  these  hydrologic  processes.  The  mechanisms  responsible  for 
sorption  include  ion  exchange,  cation  bridging,  charge  transfer,  H+ 
bonding  and  van  der  Waals  interactions  (Kookana  et  al.  1998). 

In  NSM,  physical  transport  of  inorganic  nutrients  (N  and  P)  in  surface 
runoff  is  subdivided  into  dissolved  nutrients  and  sediment-associated 
nutrients.  Thus,  there  are  two  unknown  concentrations  for  each  species  in 
surface  runoff;  the  two  unknowns  require  two  equations,  a  mass  balance 
equation  for  the  water  column  and  another  for  the  soil  layer.  For 
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suspended  solids,  a  mass  balance  equation  for  the  water  column  is 
required. 

Mathematically  the  total  concentration  of  a  species  in  the  water  column  C 
[M/L.3]  is  separated  into  two  components: 

C  =  Cd+Cp  (1) 


where: 

Cd  =  dissolved  concentration  of  a  species  [M/L3] 

Cp  =  particulate  concentration  of  a  species  [M/L3]. 

Given  the  total  concentration  and  the  two  phase  fractions,  the  dissolved 
and  particulate  concentrations  in  the  water  column  can  be  determined  as 
follows: 


Q  =  fdc 

(2) 

N  N 

c„  =  £c;  =  £/;c 

(3) 

n= 1  n= 1 


where  fd,  fp  are  fraction  of  the  total  concentration  that  are  in  dissolved 

form  and  particulate  form  associated  with  particle  “n”,  respectively.  The 
sum  of  the  two  fractions  (Cd  and  Cp)  is  equal  to  l,  and  Cp  is  the  particulate 

concentration  of  a  species  associated  with  particle  “n”  [M/L3], 

The  fraction  dissolved  and  fraction  particulate  in  the  water  column  are 
derived  using  the  distribution  coefficients  for  reversible,  linear,  equilib¬ 
rium  partitioning  between  dissolved  and  sorbed  phases  and  are  shown  in 
Equations  4  and  5  below: 


fd  = 


N 


i+Y,k„c: 


n= 1 


V  C11 
fn  _ _  -f'-d'-'s 

J  D 


N 


1+Y.K.f: 


n= 1 


(4) 


(5) 
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with  Kd  [L.3/M]  being  the  mass  transfer  coefficient.  C"  is  the  concentra¬ 
tion  of  suspended  particle  “n”  in  surface  water  [M/L3], 

For  the  soils  or  bed  sediments,  the  fractions  associated  with  dissolved  and 
particulate,  respectively,  are  derived  by  using  porosity: 

f d 2  <P  +  Kd2(X~<p)ps  (6' 

f  =  Kd2(l-<p)ps 

p2  cp+is:d2(i-cp)ps 

where: 

fd2  =  fraction  of  the  constituent  dissolved  in  the  soil  layer  (sediment 
bed) 

cp  =  porosity  of  the  upper  soil  (bed  sediment)  layer 
Kd,  Kdi  =  mass  distribution  coefficients  for  the  constituent  in  the  water 
column  and  soil  layer,  respectively  [L.3/M] 
ps  =  dry  soil  (sediment)  density  [M/L.3] 
fP2  =  fraction  of  the  constituent  in  particulate  form  in  the  soil  layer 
(sediment  bed). 
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3  Mass  Transfer  of  Dissolved  Nutrients 

Mass  transfer  flux 

Within  the  soil  or  sediment  domain,  diffusion  processes  may  dominate, 
while  at  the  soil  or  sediment  surface,  advection  processes  are  expected  to 
be  the  dominant  processes.  The  complicated  nature  of  the  flux  at  an 
interface  such  as  this  is  usually  characterized  through  the  use  of  a  mass 
transfer  coefficient,  an  empirical  coefficient  that  relates  the  concentration 
gradient  to  mass  transport  (Choy  and  Reible  1999). 

The  transfer  rate  of  dissolved  species  from  the  soil  to  the  water  column  is 
affected  by  the  concentration  gradient  across  the  water-soil  interface  as 
well  as  flow  conditions  in  the  water  column.  Mass  transfer  theory  states 
that  the  mass  flux  of  a  given  species  under  a  given  set  of  flow  conditions 
can  be  expressed  as: 


Sd  —  K{Cd2/(p  cd) 


(8) 


where: 

Sd  =  mass  transfer  flux  of  a  dissolved  species  [ML^T1] 
ke  =  mass  transfer  coefficient  between  water  column  and  sediment 
depositions  [L/T] 

Cd2  =  dissolved  concentration  of  a  species  in  the  soil  layer  [M/L3], 

The  dissolved  concentration  in  the  deposits  is  described  in  terms  of  mass 
of  the  substance  per  bulk  volume  of  the  soil  layer.  In  order  to  calculate  the 
concentrations  of  water  quality  variables  in  the  pore  water,  the  concentra¬ 
tions  of  the  substance  of  the  soil  layer  must  be  divided  by  the  porosity  as 
indicated  in  Equation  8. 

Mass  transfer  coefficient 

Overland  flow 

Surface-applied  or  soil-incorporated  chemicals  are  often  transferred  in 
significant  quantities  from  the  soil  to  surface  runoff  as  a  result  of  rainfall- 
runoff  processes.  Several  dynamic,  inter-related  processes  control 
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chemical  transport  from  the  soil  to  runoff  during  rainfall,  including  ejec¬ 
tion  of  soil  water  by  raindrop  impacts,  diffusion  of  chemicals  from  pore 
water,  infiltration,  and  soil  water-runoff  mixing  (Zhang  et  al.  1997;  Gao  et 
al.  2004).  From  early  experiments  and  calculations,  it  has  been  concluded 
that  only  a  certain  thin  soil  layer  interacts  with  the  rainfall  and  overland 
flow.  Effective  soil  depth  of  interaction  is  related  to  the  degree  of  soil 
aggregation  and  it  increases  with  land  slope,  kinetic  energy  of  raindrops, 
and  rainfall  intensity  (Ahuja  and  Lehman  1983;  Snyder  and  Woolhiser 

1985)- 


Gao  et  al.  (2004)  developed  a  model  that  combined  the  chemical  transfer 
associated  with  the  raindrop  impacts  and  diffusion  by  assuming  raindrop 
and  diffusion  processes  could  be  coupled  by  superposition.  This  model 
captured  soil-runoff  chemical  transfer  behavior  more  realistically  than 
either  mixing-layer  models  or  diffusion-based  models.  From  this  model, 
the  effective  mass  transfer  coefficient  can  be  expressed  as: 


k 


e 


(9) 


where: 


k 


m 


a 

i 

0 

Pb=  (1  — f)Ps 


diffusive  mass  transfer  coefficient  [L/T],  which  was 

derived  by  the  concentration  gradient  across  the 

hydrodynamic  boundary  layer  (Wallach  et  al.  1988, 1989) 

soil  detachability  [M/L3] 

rainfall  intensity  rate  [L/T] 

volumetric  water  content 

soil  bulk  density  [M/L3], 


Channel  flow 

Bed  release  is  an  important  source  of  chemical  diffusion  from  bed  sedi¬ 
ment  to  the  water  column.  The  mass  transfer  across  the  sediment- water 
interface  is  actually  a  diffusive  flow  that  is  dependent  upon  the  concen¬ 
tration  gradient  of  the  dissolved  mass  between  the  water  column  and  the 
benthic  sediment.  The  mass  transfer  coefficient  above  bed  sediments  was 
summarized  as  follows  by  Choy  and  Reible  (1999): 
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ke  =  0.664Re0,5  5C1/3  — ,  for  laminar  flow  boundary  layer  (10) 

ke  =  0.036  Re0  8  Sc1/3  — ,  for  turbulent  flow  boundary  layer  (11) 

Lv 


where: 

Re  =  Reynolds  number 

Sc  =  Schmidt  number 

De  =  effective  diffusion  coefficient  [L2/T] 

Lx  =  characteristic  length  across  the  sediment- water  interface  [L]. 


ERDC/EL  TR-08-25 


10 


4  Erosion  and  Sedimentation  of  Particulate 
Nutrients 

The  erosion  and  sedimentation  of  sediments  and  associated  pollutants 
shown  schematically  in  Figure  l  are  two  important  processes  in  water 
quality  modeling.  Knowledge  of  erosion  and  sedimentation  is  essential  to 
understanding  the  transport  of  pollutants  in  surface  waters.  Sediment 
detachment  by  surface  runoff  is  usually  simulated  in  terms  of  a  general¬ 
ized  erosion-deposition  theory  proposed  by  Smith  et  al.  (1995).  This 
assumes  that  the  transport  capacity  concentration  of  the  runoff  reflects  a 
balance  between  the  two  continuous  counteracting  processes  of  erosion 
and  deposition.  Some  of  these  processes  are  beyond  the  scope  of  this 
discussion;  therefore,  emphasis  will  be  placed  on  erosion  and  sedimenta¬ 
tion  as  it  relates  to  nutrients  (N  and  P).  In  general,  the  insoluble  forms  of 
nitrogen  and  phosphorus  far  exceed  their  soluble  forms.  The  physical 
transport  rate  of  both  inorganic  and  organic  forms  of  nitrogen  and 
phosphorus  with  sediment  is  described  as  follows. 


Figure  1.  Soil  layer  (sediment  bed)  and  water  interactions  for  solids. 
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Erosion 

Soil  erosion  involves  detachment,  transport,  and  subsequent  deposition. 
Sediment,  defined  here  as  eroded  soil,  is  transported  down  gradient  pri¬ 
marily  by  flowing  water.  Absorbed  on  the  surface  of  the  sediment  particles 
are  chemicals,  nutrients,  pesticides,  and  heavy  metals.  These  constituents 
are  eventually  deposited  into  nearby  streams  and  rivers  with  the  entire 
process  being  referred  to  as  sedimentation.  The  products  of  erosion  may 
be  deposited  immediately  below  their  sources,  or  may  be  transported 
considerable  distances  to  be  deposited  in  channels,  floodplains,  or  other 
water  bodies.  In  general,  the  coarser  sediments  settle  first  and  the  finer 
particles  move  further.  As  a  consequence,  suspended  sediments  tend  to 
contain  more  fines  and  less  dense  particles.  Because  these  fine  particles 
usually  contain  higher  amounts  of  adsorbed  chemicals,  it  follows  that  an 
enrichment  of  any  sorbed  chemical  in  eroded  sediment  must  be  con¬ 
sidered  and  included  in  water  quality  modeling.  Erosion  of  the  land  sur¬ 
face  is  the  major  source  of  solid  phase  pollutants  in  surface  waters.  Within 
a  watershed,  the  erosion  component  is  divided  into  upland  (hillslope) 
erosion  and  channel  erosion. 

Upland  erosion 

Soil  erosion  maybe  considered  a  natural  geomorphic  process.  In  perform¬ 
ing  an  analysis,  there  are  a  number  of  erosion  mechanisms  that  need  to  be 
taken  into  account.  There  are  types  of  erosion  in  which  flowing  water  and 
water  as  raindrops  are  directly  involved  in  the  detachment  and  movement 
of  sediments  as  well  as  a  broad  class  of  erosion  in  which  water  plays  only 
an  indirect  role  of  weakening  the  soil  matrix. 

A  common  description  of  the  soil  erosion  process  involves  detachment  of 
soil  particles  by  raindrop  impact,  transport  of  this  material  to  rills  by  sheet 
flow,  and  further  detachment  and  transport  of  soil  particles  by  the  rill  flow 
(e.g.,  Foster  et  al.  1977).  The  transport  of  sediment  by  overland  flow  is  self¬ 
regulating,  i.e.,  soil  particle  detachment  by  overland  flow  does  not  occur 
unless  excess  energy  is  available  in  addition  to  the  amount  required  to 
transport  suspended  sediments.  However,  detachment  by  rainfall  impact 
often  occurs  when  there  is  little  or  no  flow  available  for  transport.  Major 
factors  that  affect  raindrop  detachment  are  rainfall  characteristics,  soil 
characteristics,  ground  and  canopy  cover,  and  surface  water  depth 
(Taskinen  and  Bruen  2007). 
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Channel  erosion 

The  removal  of  soil  from  the  banks  or  beds  of  streams  results  in  stream 
channel  erosion.  Stream  channels  erode  as  water  flows  over  the  side  of  the 
stream  bank  or  scours  below  the  water  surface,  especially  during  severe 
floods.  The  method  for  estimating  quantitative  channel  erosion  is  the  same 
as  that  used  for  upland  erosion.  Main  differences  are  that  soil  detachment 
by  raindrop  impact  within  the  channel  is  neglected  and  that  lateral  inflows 
of  sediment  from  the  upland  become  important.  Sediment  routing  within 
channels  involves,  as  for  overland  flow,  a  comparison  of  the  availability  of 
sediment  with  the  transport  capacity  of  the  flow. 

Erosion  rate  of  particulate  nutrients 

Formulations  within  NSM  have  been  structured  such  that  the  water 
engines  provide  the  required  flow  and  sediment  fluxes  necessary  to  drive 
the  fate  and  transport  of  the  nutrients. 

The  transport  of  nutrient  particulates  from  the  soil  surface  to  the  water 
column  via  erosion  occurs  at  a  rate  that  is  proportional  to  the  rate  at  which 
sediment  particles  are  eroded  (resuspended). 

sr=E/;«r"c2  (i2> 

n= 1 


where: 

Sr  =  total  erosion  rate  of  a  nutrient  [ML  ^T1] 
v’!  =  erosion  velocity  associated  with  particle  “n”  [L/T] 

C’S2  =  particulate  concentration  associated  with  particle  “n”  in  the 
soil  layer  [M/L3], 

Sediments  in  surface  water  and  deposition 

Sediment  characteristics 

Sediments  are  usually  associated  with  streams,  rivers,  bays,  bayous,  and 
beaches.  Although  it  would  seem  that  soil  and  sediment  are  the  same,  they 
tend  to  be  very  different  in  their  physical  and  chemical  makeup,  causing 
their  interaction  with  contaminates  to  be  different.  Sediments  are  less 
variable  in  particle  size  and  often  contain  more  organic  matter  than  soil. 
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In  addition,  sediments  in  surface  water  are  contaminated  with  chemicals 
that  sorb  onto  fine-grained  organic  and  inorganic  soil  particles. 

Sediment  is  composed  of  many  materials,  including  individual  primary 
particles,  aggregates,  organic  materials,  and  associated  chemicals.  Par¬ 
ticles  can  be  mineral  or  organic  in  origin.  Mineral  particles  are  derived 
from  geologic  materials.  The  density  of  many  mineral  particles  ranges 
between  2  and  3  g/cm3  and  is  often  approximated  as  2.65  g/cm3.  Organic 
particles  are  derived  from  plant  materials,  dead  bacterial  or  algal  cells,  and 
decaying  aquatic  organisms.  Organic  particles  usually  have  a  density  only 
slightly  greater  than  that  of  water  and  contain  a  high  fraction  of  organic 
carbon,  an  excellent  sorbent  for  many  pollutants. 

The  role  of  sediment  in  water  quality  pollution  is  tied  both  to  the  particle 
size  of  sediment,  and  to  the  amount  of  particulate  organic  carbon  associ¬ 
ated  with  the  sediment.  Size,  shape,  and  density  affect  the  settling  velocity, 
which  in  turn  affects  sediment  transport  rates  and  sediment  deposition. 
For  phosphorus  and  metals,  particle  size  (usually  63  pm  (silt  +  clay)  frac¬ 
tion)  is  of  primary  importance  due  to  the  large  surface  area  of  very  small 
particles.  Phosphorus  and  metals  tend  to  be  highly  attracted  to  ionic 
exchange  sites  that  are  associated  with  clay  particles  and  with  the  iron  and 
manganese  coatings  that  commonly  occur  on  these  small  particles.  Many 
of  the  persistent,  bioaccumulating  and  toxic  organic  contaminants,  espe¬ 
cially  chlorinated  compounds  including  many  pesticides,  are  strongly 
associated  with  sediment  and  especially  with  the  organic  carbon  that  is 
transported  as  part  of  the  sediment  load  in  surface  water. 

Deposition  rate  of  particulate  nutrients 

The  magnitude  of  the  deposition  flux  of  a  contaminant  is  equal  to  the 
product  of  the  rate  of  sediment  deposition  and  the  contaminant  concen¬ 
tration  associated  with  the  settling  particles.  Settling  velocity  depends  not 
only  on  the  size,  shape,  and  density  of  particles,  but  also  on  the  concen¬ 
tration  of  the  particles  provided  by  the  water  engines. 

The  deposition  of  particulate  nutrients  from  the  water  column  is  expressed 
as: 


n 

s 


n= 1 


(13) 
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where: 

Ss  =  total  deposition  rate  of  a  nutrient  [ML^T1] 

Vg  =  settling  velocity  associated  with  particle  “n”  [L/T] 

Cg  =  particulate  concentration  associated  with  particle  “n”  in  the 
water  column  [M/L3], 
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5  Soil  Plant  Dynamics 

In  order  to  accurately  simulate  the  removal  of  water  and  nutrients  from 
the  plant  root  zone,  transpiration,  and  biomass/yield  production,  a 
dynamic  soil  plant  module  has  been  developed  for  NSM.  The  soil  plant 
algorithms  have  been  taken  from  the  SWAT/EPIC  model  (Neitsch  et  al. 
2002).  Since  its  inception,  EPIC  has  evolved  into  a  comprehensive  agro¬ 
ecosystem  model  capable  of  simulating  the  growth  of  crops  in  complex 
rotations  and  management  operations,  such  as  tillage,  irrigation,  fertiliza¬ 
tion  and  liming  (Williams  et  al.  1984).  This  chapter  provides  a  detailed 
description  of  the  plant  growth  processes  that  can  be  found  within  the 
modules. 

Currently  NSM  can  simulate  approximately  too  plant  species  including 
crops,  native  grasses,  and  trees.  The  process  descriptions  utilize  a  single 
plant  growth  model  to  simulate  all  types  of  land  covers.  The  algorithms  are 
able  to  differentiate  between  annual  and  perennial  plants.  Annual  plants 
grow  from  the  planting  date  to  the  harvest  date  or  until  the  accumulated 
heat  units  equal  the  potential  heat  units  for  the  plant.  Perennial  plants 
maintain  their  root  systems  throughout  the  year,  becoming  dormant  in  the 
winter  months  and  resume  growth  when  the  average  daily  air  temperature 
exceeds  the  minimum,  or  base,  temperature  required. 

For  each  increment  of  time,  potential  plant  growth,  i.e.  plant  growth  under 
ideal  growing  conditions,  is  calculated.  Ideal  growing  conditions  consist  of 
adequate  water  and  nutrient  supply  and  a  favorable  climate.  Differences  in 
growth  between  plant  species  are  defined  by  the  parameters  contained  in 
the  plant  growth  database  (Neitsch  et  al.  2002). 

Plant  growth 

Plant  growth  is  modeled  by  simulating  leaf  area  development,  light  inter¬ 
ception,  and  conversion  of  intercepted  light  into  biomass  assuming  a  plant 
species-specific  radiation-use  efficiency  (RUE). 
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Potential  growth 

Biomass  production 

NSM  uses  the  concept  of  RUE  by  which  a  fraction  of  daily  photosyn- 
thetically  active  radiation  (PAR)  is  intercepted  by  the  plant  canopy  and 
converted  into  plant  biomass.  Daily  gains  in  plant  biomass  are  affected  by 
air  vapor  pressure  deficits  and  atmospheric  CO2  concentration  (Stockle  et 
al.  1992).  Stress  indices  for  water,  temperature,  N,  P,  and  aeration  are 
calculated  daily  using  the  value  of  the  most  severe  stress  to  reduce  poten¬ 
tial  plant  growth  and  crop  yield.  Similarly,  stress  factors  for  soil  strength, 
temperature,  and  aluminum  toxicity  are  used  to  adjust  potential  root 
growth  (Jones  et  al.  1991). 


The  amount  of  daily  solar  radiation  intercepted  by  the  leaf  area  of  the 
plant  is  calculated  using  the  Beer’s  Law  (Chapra  1997): 


R 


phosyn 


0-5  Rday 


1-exp  (krLAl) 


(14) 


where: 

Rphosyn  =  amount  of  intercepted  PAR  on  a  given  day  (M  J  nr2) 

Rday  =  incident  total  solar  radiation  (MJ  nr2) 

0.5 Rday  =  incident  PAR  (MJ  nr2) 

ki  =  light  extinction  coefficient,  which  is  taken  as  -0.65  for  all 
plants 

LAI  =  leaf  area  index  for  a  given  day. 


PAR  is  radiation  with  a  wavelength  between  400  and  700  mm  (McCree 
1972).  Direct  solar  beam  radiation  contains  roughly  45  percent  of  PAR 
while  diffuse  radiation  contains  approximately  60  percent  of  PAR 
(Monteith  1972).  The  fraction  of  PAR  will  vary  from  day  to  day  with 
variation  in  cloud  cover. 


RUE,  the  amount  of  dry  biomass  produced  per  unit  intercepted  solar  radi¬ 
ation,  is  defined  in  the  plant  growth  database,  assumed  to  be  independent 
of  the  plant’s  growth  stage.  The  maximum  increase  in  biomass  on  a  given 
day  that  will  result  from  the  intercepted  PAR  is  assumed  to  be  a  linear 
function  of  radiant  energy  (Monteith  1977): 
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Abio  =  RUE  ■  Rphosyn 


(15) 


where: 

A  bio  =  potential  increase  in  total  plant  biomass  on  a  given  day 
(kg/ha) 

RUE  =  RUE  of  the  plant  (kg/hg  (MJ/m2)'1  or  to1  g/MJ). 

The  total  biomass  on  a  given  day  d  is  calculated  as: 

d 

bio  =  y^Abipj  (16) 

i=i 

where: 

bio  =  total  plant  biomass  on  a  given  day  (kg/ha) 

A bioi  =  increase  in  total  plant  biomass  on  a  day,  i  (kg/ha). 

Impact  of  climate  on  radiation-use  efficiency.  RUE  is  sensitive  to 
variations  in  atmospheric  CO2  concentrations  and  equations  have  been 
incorporated  into  the  EPIC  formulations  to  modify  the  default  RUE  values 
in  the  plant  database  for  climate  change  studies.  The  relationship  used  to 
adjust  the  RUE  for  effects  of  elevated  CO2  is  given  by  Stockle  et  al.  (1992): 


RUE 


100  C02 

C02  +  exp  (r±  -  r2C02 ) 


(17) 


where: 


C02 

ri  and  r2 
ri 
T2 


concentration  of  carbon  dioxide  in  the  atmosphere  (parts  per 

million  by  volume,  or  ppmv) 

shape  coefficients 

first  shape  coefficient 

second  shape  coefficient. 


The  shape  coefficients  are  calculated  by  solving  Equation  17  using  two 
known  points  ( [RUEamb ,  C02amb)  and  (RUE hi,  C02hi ): 
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r±  =  In 


CO 


2amb 


0.01  RUEamb) 


CO. 


2  amb 


+  pCO 


2  amb 


(18) 


In 


CO. 


2  amb 


0.01RUEamb) 


CO 


2  amb 


In 


CO. 


2  hi 


0.01  RUEhi) 


CO 


2  hi 


C02hi  C02amh 


(19) 


where: 


C02amb 

RUEamb 


COlhi 

RUEhi 


ambient  atmospheric  CO2  concentration  (ppmv) 

RUE  of  the  plant  at  ambient  atmospheric  CO2  concentration 

(kg/hg  (MJ/m2)-1  or  ter1  g/MJ) 

elevated  atmospheric  CO2  concentration  (ppmv) 

RUE  of  the  plant  at  the  elevated  atmospheric  CO2 
concentration,  CO2M,  (kg/hg  (MJ/m2)-1  or  to1  g/MJ). 


Equation  17  was  developed  when  the  ambient  atmospheric  CO2  concentra¬ 
tion  was  330  ppmv  and  is  valid  for  CO2  concentrations  in  the  range  of 
330-660  ppmv.  Even  though  the  ambient  atmospheric  concentration  of 
C02  is  now  higher  than  330  ppmv,  this  value  is  still  used  in  the  calculation. 
If  the  CO2  concentration  used  in  the  simulation  is  less  than  330  ppmv,  the 
model  defines  RUE  =  RUEamb. 


Stockle  and  Kiniry  (1990)  have  shown  that  a  plant’s  RUE  is  affected  by  air 
vapor  pressure  deficit.  For  a  plant,  a  threshold  vapor  pressure  deficit  is 
defined  at  which  the  plant’s  RUE  begins  to  drop  in  response  to  the  vapor 
pressure  deficit.  The  adjusted  RUE  is  calculated: 

RUE  =  RUEvpd=1  -  A ruedcl  ( vpd  -  vpdthr )  if  vpd  >  vpdthr  (20) 
RUE  =  RUEvpd=1  if  vpd  <  vpdthr  (21) 


where: 


RUEvpd= 1  =  RUE  for  the  plant  at  a  vapor  pressure  deficit  of  1  kPa  (kg/hg 
(MJ/m2)-1  or  to1  g/MJ) 

A ruedd  =  rate  of  decline  in  RUE  per  unit  increase  in  vapor  pressure 
deficit  (kg/hg  (MJ/m2)1  or  to1  g/MJ) 
vpd  =  vapor  pressure  deficit  (kPa) 
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vpdthr  =  threshold  vapor  pressure  deficit  above  which  a  plant  will 
exhibit  reduced  RUE  (kPa). 

The  RUE  value  reported  for  the  plant  in  the  plant  growth  database  RUEamb 
or  adjusted  for  elevated  CO2  levels  (Equation  17)  is  the  value  used  for 
RUEvpd=  1.  The  threshold  vapor  pressure  deficit  for  reduced  RUE  is 
assumed  to  be  1.0  kPa  for  all  plants  ( vpdthr  =  1.0).  The  RUE  is  never 
allowed  to  fall  below  27  percent  of  RUEamb ■  This  minimum  value  was 
based  on  field  observations  (Neitsch  et  al.  2002). 


Canopy  cover  and  height.  In  the  initial  period  of  plant  growth,  canopy 
height  and  leaf  area  development  are  controlled  by  the  optimal  leaf  area 
development  curve: 


fii 


fii 


PHU 


LAI,  max 


fr PHU  +  eXP  (4  4  '  fi PHU  ) 


(22) 


where: 


frlAI,  max 

frpHu 
h  and  h 


fraction  of  the  plant’s  maximum  LAI  corresponding  to  a  given 

fraction  of  potential  heat  units  for  the  plant 

fraction  of  potential  heat  units  accumulated  for  the  plant  on  a 

given  day  in  the  growing  season 

shape  coefficients. 


The  fraction  of  potential  heat  units  accumulated  by  a  given  date  is 
calculated: 


EHUi 

J^PHU  -  pHU 


(23) 


where: 

HUi  =  heat  units  accumulated  on  day  i  (heat  units) 

PHU  =  total  potential  heat  units  for  the  plant  (heat  units). 

The  shape  coefficients  are  calculated  by  solving  Equation  22  using  two 
known  points  (friAi, 1,  Jtphu,  1)  and  if  run, 2, /tphu, 2 ): 
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Z±=  In 


PHU,  1 


frl 


'  frpHU,l 


LAI,  1 


+  4 'A 


PHU,  1 


Zo  = 


In 


fri 


PHU,  1 


frl 


PHU,  1 


LAI,1 


In 


PHU, 2 


frl 


frpHU,  2 


LAI, 2 


fVPHU,2  frl 


PHU,  1 


(24) 


(25) 


where: 


frpHU,  1 


Jtlai,i 


frpHU, 2 
frlAI,2 


fraction  of  the  growing  season  (i.e.,  fraction  of  total  potential 
heat  units)  corresponding  to  the  first  point  on  the  optimal  leaf 
area  development  curve 

fraction  of  the  maximum  plant  LAI  (i.e.,  fraction  of  LAImax ) 
corresponding  to  the  first  point  on  the  optimal  leaf  area 
development  curve 

fraction  of  the  growing  season  corresponding  to  the  second 
point  on  the  optimal  leaf  area  development  curve 
fraction  of  the  maximum  plant  LAI  corresponding  to  the 
second  point  on  the  optimal  leaf  area  development  curve. 


The  canopy  height,  on  a  given  day,  is  calculated: 


Zf  :  Zlc  max 


LAI,  max 


(26) 


where: 

hc  =  canopy  height  for  a  given  day  (m) 
hc, max  =  plant’s  maximum  canopy  height  (m). 

The  amount  of  canopy  cover  is  expressed  as  the  LAI.  The  leaf  area  added 
on  day  i  is  calculated: 


MAT. 


LAI,  max 


frl 


i- 1 


LAI,  max 


LAI 


max 


{l  -  exp  [5  (LAI^  -  LAImax )]} 


(27) 
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and  the  total  LAI  is  calculated: 


LAIi=LAIi_1+ALAIi 


(28) 


where: 


ALAI, 
LAI,  and  LAIi-i 

Jtlai,  max2  and  frLAi,  max2  ^ 


LAI n 


leaf  area  added  on  day  i 
leaf  area  indices  for  day  i,  and  i-i  respectively 
fraction  of  the  plant’s  maximum  LAI  calculated 
with  Equation  22  for  day  i  and  z-i 
maximum  LAI  for  the  plant. 


LAI  is  defined  as  the  area  of  green  leaf  per  unit  area  of  land  (Watson  1947). 
Once  the  maximum  leaf  area  index  is  reached,  LAI  will  remain  constant 
until  leaf  senescence  begins  to  exceed  leaf  growth.  In  general  terms,  leaf 
senescence  is  a  way  for  a  deciduous  plant  to  prepare  for  winter  and  recycle 
some  of  the  valuable  and  often  scarce  mineral  nutrients,  such  as  nitrogen 
and  phosphorus.  Leaf  senescence  is  also  a  way  to  get  rid  of  old  and  photo- 
synthetically  less  efficient  leaves  in  both  deciduous  and  evergreen  plants. 
Once  leaf  senescence  becomes  the  dominant  growth  process,  the  LAI  is 
calculated: 


LAI  =  16 LAImax  (l  -  frPHU)2  frPHu  >  frpHu,sen  (29) 

where  frpHu, sen  is  the  fraction  of  growing  season  (PHU)  at  which 
senescence  becomes  the  dominant  growth  process. 

Root  development.  The  amount  of  total  plant  biomass  partitioned  to 
the  root  system  is  30-50  percent  in  seedlings  and  decreases  to  5-20  per¬ 
cent  in  mature  plants  (Jones  1985).  Within  NSM,  the  fraction  of  total 
biomass  in  roots  varies  from  0.40  at  emergence  to  0.20  at  maturity.  The 
daily  root  biomass  fraction  is  calculated  with  the  equation: 

frroot  —  0-40  -  0. 20 frPHU  (30) 

where  fr root  is  the  fraction  of  total  biomass  partitioned  to  roots  on  a  given 
day  in  the  growing  season. 
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Calculation  of  root  depth  varies  according  to  plant  type.  NSM  assumes 
perennials  and  trees  have  roots  down  to  the  maximum  rooting  depth 
defined  for  the  soil  throughout  the  growing  season: 

7  -  7 

root  root,  max 


where: 


Zroot  =  depth  of  root  development  in  the  soil  on  a  given  day  (mm) 
Zroot,  max  =  maximum  depth  for  root  development  in  the  soil  (mm). 

The  simulated  root  depth  for  annuals  varies  linearly  from  o.o  mm  at  the 
beginning  of  the  growing  season  to  the  maximum  rooting  depth  at 
Jtphu  =  0.40  using  the  equations: 


Z root  2.5  frpHU  '  Zroot, max 

if  frpHU  <0.40 

(32) 

7  —  7 

root  roof, max 

iffr  >  0.40 

(33) 

The  maximum  rooting  depth  is  defined  by  comparing  the  maximum 
potential  rooting  depth  for  the  plant  from  the  plant  growth  database  and 
the  maximum  potential  rooting  depth  for  the  soil.  The  smaller  of  the  two  is 
taken  as  the  maximum  rooting  depth. 

Maturity.  Plant  maturity  is  reached  when  the  fraction  of  potential  heat 
units  accumulated/rpw is  equal  to  1.0.  Once  maturity  is  reached,  the  plant 
ceases  to  transpire  and  take  up  water  and  nutrients.  Simulated  plant  bio¬ 
mass  remains  stable  until  the  plant  is  harvested  or  killed  as  the  result  of  a 
management  operation. 

Actual  plant  growth 

Net  growth  occurs  in  plants  as  a  function  of  mortality,  respiration  require¬ 
ments,  photosynthesis,  and  demands  by  the  plant  organs.  The  plant 
growth  factor  quantifies  the  fraction  of  potential  growth  achieved: 

yreg  =l-max(wstrs,nstrs,pstrs)  (34) 

where  yreg  is  the  plant  growth  factor  and  ranges  from  0.0  to  1.0. 
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The  potential  biomass  predicted  is  adjusted  if  one  of  the  four  plant  stress 
factors  is  greater  than  o.o: 


Abi0act=Abi0-yreg  (35) 

The  potential  leaf  area  added  is  also  adjusted  for  plant  stress: 

ALAI  acti  =  ALAI  (36) 


Plant  growth  constraints 

The  previous  sections  discussed  potential  growth.  In  reality,  plant  growth 
may  be  constrained  due  to  extreme  temperatures,  and  insufficient  water, 
nitrogen,  or  phosphorus. 

Temperature  stress.  Temperature  stress  is  a  function  of  the  daily 
average  air  temperature  and  the  optimal  temperature  for  plant  growth.  As 
the  air  temperature  diverges  from  the  optimal,  the  plant  will  begin  to 
experience  stress.  The  equations  used  to  determine  temperature  stress  are: 

tStrS  =  1.0  f(w  <  Tbase  (37) 


tstrs 


1-exp 


■0.10541  Top,-Ta  '2 


T  —  T 

av  base  / 


\2 


Tbase  <  Tav  <  T0pt  (38) 


tstrs 


1-exp 


-o.io54(r0i„-r„)2 
1 2  T  —  T  —T  j2 

\  opt  av  base) 


T0pt  <  Tav  <  2*T0prTbase 


(39) 


tstrs  =  1.0 


Tav  >  2.*T0pt~Tbase  (40) 


where: 

tstrs  =  temperature  stress  for  a  given  day  expressed  as  a  fraction  of 
optimal  plant  growth 
Tav  =  mean  air  temperature  for  day  (°C) 

Tbase  =  plant’s  base  or  minimum  temperature  for  growth  (°C) 

Topt  =  plant’s  optimal  temperature  for  growth  (°C). 
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Water  stress.  Water  stress  is  simulated  by  comparing  actual  and 
potential  plant  transpiration: 


wstrs  =  1  -  =  1  -  WactualuP 

Et  Et 


(41) 


where: 


wstrs 

Et,act 

Et 

Wactualup 


water  stress 

actual  amount  of  transpiration  on  a  given  day  (mm  H20) 
maximum  plant  transpiration  (mm  H20) 
total  plant  water  uptake  for  the  day  (mm  H20). 


Nitrogen  stress.  Nitrogen  stress  is  calculated  only  for  non-legumes. 
Legumes  do  not  experience  nitrogen  stress  due  to  nitrogen  fixation. 


Nitrogen  stress  is  quantified  by  comparing  actual  and  optimal  plant  nitro¬ 
gen  levels.  Nitrogen  stress  varies  non-linearly  between  o.o  at  optimal 
nitrogen  content  and  l.o  when  the  nitrogen  content  of  the  plant  is 
50  percent  or  less  of  the  optimal  value.  Nitrogen  stress  is  computed  with 
the  following  equation: 


nstrs  =  1 - ' 11 - 

cpn  +exp(3.535  —  0.02597cp;!) 


(42) 


where: 

nstrs  =  nitrogen  stress  for  a  given  day 
cpn  =  scaling  factor  for  nitrogen  stress. 

The  scaling  factor  is  calculated: 


bio 


N 


bio 


N,opt 


(Pn  =  200 


0.5 


(43) 
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where: 

Mon  =  actual  mass  of  nitrogen  stored  in  plant  material  (kg  N/ha) 
bioN,opt  =  optimal  mass  of  nitrogen  stored  in  plant  material  for  the 
current  growth  stage  (kg  N/ha). 

Phosphorus  stress.  As  with  nitrogen,  phosphorus  stress  is  quantified  by 
comparing  actual  and  optimal  plant  phosphorus  levels.  Phosphorus  stress 
varies  non-linearly  between  o.o  at  optimal  phosphorus  content  and  1.0 
when  the  phosphorus  content  of  the  plant  is  50  percent  or  less  of  the  opti¬ 
mal  value.  Phosphorus  stress  is  computed  with  the  following  equation: 

(p  n 

vstrs  =  1 - — -  (44) 

cpp+ exp  (3.535 -0.02597(pp) 


where: 

pstrs  =  phosphorus  stress  for  a  given  day 
(pp  =  scaling  factor  for  phosphorus  stress. 

The  scaling  factor  is  calculated: 


<Pr 


200 


bio. 


bio 


0.5 


P,opt 


(45) 


where: 

biop  =  actual  mass  of  phosphorus  stored  in  plant  material  (kg  P/ha) 
biop,opt  =  optimal  mass  of  phosphorus  stored  in  plant  material  for  the 
current  growth  stage  (kg  P/ha). 

Water  uptake  by  plants 

A  simple  algorithm  is  used  to  simulate  multiple  layer  soil  water  flow,  plant 
water  uptake  and  plant-water  relations  for  a  given  soil-plant-atmosphere 
system.  The  schematic  representation  in  the  soil  profile  is  illustrated  in 
Figure  2. 
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The  potential  water  uptake  from  the  soil  surface  to  any  depth  in  the  root 
zone  is  estimated  with  the  equation: 


wup,z  = 


l-exp(-Pl0) 


1-exp 


-P, 


root 


(46) 


where: 

ivup,z  =  potential  water  uptake  from  the  soil  surface  to  a  specified 
depth  z  on  a  given  day  (mm  H20) 

(3U>  =  water-use  distribution  parameter 
z  =  depth  from  the  soil  surface  (mm). 

The  potential  water  uptake  from  any  soil  layer  can  be  calculated  by  solving 
Equation  46  for  the  depth  at  the  top  and  bottom  of  the  soil  layer  and 
taking  the  difference: 


^ up,ly  ^ up,zl  ^ up,zu 


(47) 
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where: 

ivUpfy  =  potential  water  uptake  for  soil  layer  ly  (mm  H20) 
wup,zi  =  potential  water  uptake  for  the  profile  at  the  lower  boundary  of 
the  soil  layer  (mm  H20) 

Wup,zu  =  potential  water  uptake  for  the  profile  at  the  upper  boundary  of 
the  soil  layer  (mm  H20). 

Since  root  density  is  greatest  near  the  soil  surface  and  decreases  with 
depth,  the  water  uptake  from  the  upper  layers  is  assumed  to  be  much 
greater  than  that  in  the  lower  layers.  The  water-use  distribution  parameter 
Pw  is  set  to  to  in  this  formulation.  With  this  value,  50  percent  of  the  water 
uptake  will  occur  in  the  upper  6  percent  of  the  root  zone. 

The  amount  of  water  uptake  that  occurs  on  a  given  day  is  a  function  of  the 
amount  of  water  required  for  the  plant  for  transpiration  Et  and  the  amount 
of  water  available  in  the  soil  SW  (mm  H20).  Equations  46  and  47  calculate 
potential  water  uptake  solely  as  a  function  of  water  demand  for  transpira¬ 
tion  and  the  depth  distribution  defined  in  Equation  46.  Within  this  formu¬ 
lation,  the  plant  dynamics  module  modifies  the  initial  potential  water 
uptake  from  a  given  soil  layer  to  reflect  soil  water  availability  in  the  follow¬ 
ing  ways. 

If  the  upper  layers  in  the  soil  profile  do  not  contain  enough  water  to  meet 
the  potential  water  uptake  calculated  in  Equation  47,  users  may  allow 
lower  layers  to  compensate.  The  equation  used  to  calculate  the  adjusted 
potential  water  uptake  is: 

Wup, ly  =  Wup,ly  +  W demand  '  P™  (48) 


where: 

w'upjy  =  adjusted  potential  water  uptake  for  layer  ly  (mm  H20) 

Wdemand  =  water  uptake  demand  not  met  by  overlying  soil  layers  (mm 
H20) 

epco  =  plant  uptake  compensation  factor. 


The  plant  uptake  compensation  factor  can  range  from  o.oi  to  1.00  and  is 
set  by  the  user.  As  epco  approaches  l.o,  the  model  allows  more  of  the 
water  uptake  demand  to  be  met  by  lower  layers  in  the  soil.  As  epco 
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approaches  o.o,  the  model  allows  less  variation  from  the  depth  distri¬ 
bution,  described  in  Equation  46,  to  take  place. 

As  the  water  content  of  the  soil  decreases,  the  water  in  the  soil  is  held 
more  and  more  tightly  by  the  soil  particles  and  it  becomes  increasingly 
difficult  for  the  plant  to  extract  water  from  the  soil.  To  reflect  the  decrease 
in  the  efficiency  of  the  plant  in  extracting  water  from  dryer  soils,  the 
potential  water  uptake  is  modified  using  the  following  equations: 


w 


miy 


=  w 


upjy 


exp 


SW, 


iy 


0.25  AWC 


iy 


SWiy<  (0.25AWCiy)  (49) 


^ up,ly  ^ up,ly 


SWiy  >  (0.25 AWCly)  (50) 


where: 

w  'm  ly  =  potential  water  uptake  adjusted  for  initial  soil  water  content 
(mm  H20) 

SWiy  =  amount  of  water  in  the  soil  layer  on  a  given  day  (mm  H20) 
AWCiy  =  available  water  capacity  for  layer  ly  (mm  H20). 


The  available  water  capacity  is  calculated: 


AWCly  =  FCly-WPly 


(51) 


where: 

FCiy  =  water  content  of  layer  ly  at  field  capacity  (mm  H20) 
WPiy  =  water  content  of  layer  ly  at  wilting  point  (mm  H20). 


Once  the  potential  water  uptake  has  been  modified  for  soil  water  condi¬ 
tions,  the  actual  amount  of  water  uptake  from  the  soil  layer  is  calculated: 


w 


actualup,  ly 


=  mm 


w 


miy’ 


SWly-WPly 


(52) 


where  iv actualup, iy  is  the  actual  water  uptake  for  layer  ly  (mm  H20). 
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The  total  water  uptake  for  the  day  is  calculated: 


w 


actualup 


=z> 

ly= i 


actualup,  ly 


(53) 


where  n  is  the  number  of  layers  in  the  soil  profile. 

The  total  plant  water  uptake  for  the  day,  calculated  with  Equation  53,  is 
also  the  actual  amount  of  transpiration  that  occurs  on  the  day. 

^t,  act  =  W actualup  (54) 


Nutrient  uptake  by  plants 

The  current  formulations  account  for  nitrogen  and  phosphorus  uptake  by 
plants.  The  simulation  of  nitrogen  and  phosphorus  uptake  and  partition¬ 
ing  to  plants  follows  the  approach  of  the  EPIC  model.  The  approach  allows 
simulation  of  changing  allocations  of  nitrogen  and  phosphorus  based  on 
availability  and  demand. 


Nitrogen  uptake 

Nitrogen  is  taken  into  plant  roots  from  soil  solution  as  nitrate  (IV03). 
Ammonium  ( NH4 )  can  also  be  used  as  a  source  of  nitrogen,  but  must  first 
be  converted  to  nitrate  by  soil  micro-organisms.  Plant  nitrogen  uptake  is 
controlled  by  the  plant  nitrogen  equation.  The  plant  nitrogen  equation 
calculates  the  fraction  of  nitrogen  in  the  plant  biomass  as  a  function  of 
growth  stage  given  optimal  growing  conditions: 


frN  —  {frN,  i  frN, : 


fri 


PHU 


frpmj  +  exp  ( n±  n2frPHU ) 


+  fri 


N,3 


(55) 


where: 


Jtn  = 

frN,  1  = 
frN, 3  = 

ni  and  ri2= 


fraction  of  nitrogen  in  the  plant  biomass  on  a  given  day 
normal  fraction  of  nitrogen  in  the  plant  biomass  at  emergence 
normal  fraction  of  nitrogen  in  the  plant  biomass  at  maturity 
shape  coefficients  calculated  by  solving  Equation  55  using  two 
known  points  (fr n, 2,  frpHu, 50%)  and  (frjv,3,/rpw,ioo%): 


ERDC/EL  TR-08-25 


30 


n±  =  In 


fi'l 


PHU,  50% 


fi'N,  2  (fi’l 


N,  3 


■fi'l 


PHU,  50% 


fi'N,  i 


+  «2>J 


PHU,  50% 


(56) 


n,  = 


In 


PHU,  50% 


1- 


frN,2~fil 

-fi'l 


N,3 


~  fri 


PHU,  50% 


N,3 


-In 


fi'l 


PHU,  100% 


1- 


fi'N,~3~fi'i 


N,3 


-fi'l 


PHU,  100% 


{fi'N,!  fi'N, 3) 


fi'l 


PHU,  100% 


PHU, 50% 


(57) 


where: 


fr phu, 50%  =  fraction  of  potential  heat  units  accumulated  for  the  plant  at 
50  percent  maturity  (fr phu, 50%  =  0.50) 
fr phu,  100%=  fraction  of  potential  heat  units  accumulated  for  the  plant  at 
maturity  (frpHu,  100%  =  1.0). 


The  normal  fraction  of  nitrogen  in  the  plant  biomass  near  maturity  is  used 


in  Equation  57  to  ensure  that  the  denominator  term 
does  not  equal  1.  The  formulation  assumes  (frN>„3  - 


^  (frN,~3  J^N,  3) 
( fi'N, ,1  —  fi'N,  3) 

frN,  3)  =  0.00001. 


To  determine  the  mass  of  nitrogen  that  should  be  stored  in  the  plant 
biomass  on  a  given  day  the  nitrogen  fraction  is  multiplied  by  the  total 
plant  biomass: 


bioN,oPt=fi'N-bio  (58) 

The  plant  nitrogen  demand  for  a  given  day  is  determined  by  taking  the 
difference  between  the  nitrogen  content  of  the  plant  biomass  expected  for 
the  plant’s  growth  stage  and  the  actual  nitrogen  content: 

Nup  =  bio  N,  opt  —  bi°N  (59) 


where  Nup  is  the  potential  nitrogen  uptake  (kg  N/ha). 
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The  depth  distribution  of  nitrogen  uptake  is  calculated  with  the  equation: 


N, 


up 


N ri-exp(-P„)l 


1-exp 


-p„- 


root 


(60) 


where: 

Nup,z  =  potential  nitrogen  uptake  from  the  soil  surface  to  depth  z 
(kg  N/ha) 

(3n  =  nitrogen  uptake  distribution  parameter. 

The  potential  nitrogen  uptake  for  a  soil  layer  is  calculated  by  solving 
Equation  60  for  the  depth  at  the  upper  and  lower  boundaries  of  the  soil 
layer  and  taking  the  difference: 

N  up,  ly  =  ^  up,  zl  —  N up,  ZU  (61) 


where: 

Nup,iy  =  potential  nitrogen  uptake  for  layer  ly  (kg  N/ha) 

Nup,zi  =  potential  nitrogen  uptake  from  the  soil  surface  to  the  lower 
boundary  of  the  soil  layer  (kg  N/ha) 

Nup,ZU  =  potential  nitrogen  uptake  from  the  soil  surface  to  the  upper 
boundary  of  the  soil  layer  (kg  N/ha). 

Since  root  density  is  greatest  near  the  surface,  nitrogen  uptake  in  the 
upper  portion  of  the  soil  will  be  greater  than  in  the  lower  portion.  The 
depth  distribution  of  nitrogen  uptake  is  controlled  by  pn,  the  nitrogen 
uptake  distribution  parameter,  a  variable  users  are  allowed  to  adjust. 

Nitrogen  removed  from  the  soil  by  plants  is  taken  from  the  nitrate  pool. 
The  importance  of  the  nitrogen  uptake  distribution  parameter  lies  in  its 
control  over  the  maximum  amount  of  nitrate  removed  from  the  upper 
layers.  Because  the  top  10  mm  of  the  soil  profile  interacts  with  surface  run¬ 
off,  the  nitrogen  uptake  distribution  parameter  will  influence  the  amount 
of  nitrate  available  for  transport  in  surface  runoff.  This  formulation  allows 
lower  layers  in  the  root  zone  to  fully  compensate  for  lack  of  nitrate  in  the 
upper  layers,  so  there  should  not  be  significant  changes  in  nitrogen  stress 
with  variation  in  the  value  used  for  (3„. 
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The  actual  amount  of  nitrogen  removed  from  a  soil  layer  is  calculated: 

Nactualup,ly  =  up, ly  +  N  demand ’N03ly)  (62) 

where: 

NactuaiuP,iy=  actual  nitrogen  uptake  for  layer  ly  (kg  N/ha) 

Ndemand  =  nitrogen  uptake  demand  not  met  by  overlying  soil  layers 
(kg  N/ha) 

N03iy  =  nitrate  content  of  soil  layer  ly  (kg  N03-N/ha). 

Nitrogen  fixation 

If  nitrate  levels  in  the  root  zone  are  insufficient  to  meet  the  demand  of  a 
legume,  this  formulation  allows  the  plant  to  obtain  additional  nitrogen 
through  nitrogen  fixation.  Nitrogen  fixation  is  calculated  as  a  function  of 
soil  water,  soil  nitrate  content,  and  growth  stage  of  the  plant: 

Nfix  ^ demand  f gr  '  (  f sw  ’  f N03  ’  1)  (63) 

where: 


Nnx  =  amount  of  nitrogen  added  to  the  plant  biomass  by  fixation 
(kg  N/ha) 

fgr  =  growth  stage  factor  (o.o  -  l.o) 
fsw  =  soil  water  factor  (o.o  -  l.o) 
fm3  =  soil  nitrate  factor  (o.o  -  l.o). 

The  maximum  amount  of  nitrogen  that  can  be  fixed  by  the  plant  on  a  given 

day  is  Ndemand- 

Growth  stage  exerts  the  greatest  impact  on  the  ability  of  the  plant  to  fix 
nitrogen.  The  growth  stage  factor  is  calculated: 


fgr 

U 

fgr 


0.0 

frpHu  <0.15 

(64) 

6-67  frpHu -1-0 

0.15  <  frpHu  <  0.30 

(65) 

1.0 

0.30  <  frpHu  ^  0.55 

(66) 
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fgr  —  3.75  5.0  frpHU 

0.55  <  Jtphu  ^  0.75 

(67) 

O 

o 

SL, 

frpHU  >  0.75 

(68) 

The  soil  nitrate  factor  inhibits  nitrogen  fixation  as  the  presence  of  nitrate 
in  the  soil  goes  up.  The  soil  nitrate  factor  is  calculated: 

fN0 3=1-0  N03  <  100  (69) 

fN0 3  =  1.5  —  0.00  05AT03  100  <  N03  <  300  (70) 

fm 3=0.0  N03  >  300  (71) 

where  N03  is  the  nitrate  content  of  the  soil  profile  (kg  N03-N/ha). 

The  soil  water  factor  is  calculated: 


SW 

O.850/c 


(72) 


where: 

SIT  =  amount  of  water  in  the  soil  profile  (mm  H20) 

0/c  =  water  content  of  the  soil  profile  at  field  capacity  (mm  H20). 

Phosphorus  uptake 

Phosphorus  is  an  important  plant  macronutrient,  making  up  about 
0.2  percent  of  a  plant’s  dry  weight.  It  is  a  component  of  key  molecules  such 
as  nucleic  acids,  phospholipids,  and  ATP,  and,  consequently,  plants  cannot 
grow  without  a  reliable  supply  of  this  nutrient.  Plant  phosphorus  uptake  is 
controlled  by  the  plant  phosphorus  equation.  The  plant  phosphorus  equat¬ 
ion  calculates  the  fraction  of  phosphorus  in  the  plant  biomass  as  a  function 
of  growth  stage  given  optimal  growing  conditions: 


frp  ~  ifrp.i  ~  fri 


'  P,  3 


frl 


PHU 


frPHU  +  exp(  Pi  p2  *  fr PHU  ) 


+/0 


P,  3 


(73) 
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where: 


frp  =  fraction  of  phosphorus  in  the  plant  biomass  on  a  given  day 
frp, i  =  normal  fraction  of  phosphorus  in  the  plant  biomass  at 


emergence 

frp, 3  =  normal  fraction  of  phosphorus  in  the  plant  biomass  at 
maturity 

pi  and  p2  =  shape  coefficients  calculated  by  solving  Equation  73  using  two 
known  points  {frp, 2,  frpHu, 50%)  and  (frp, 3,  Jtphu,  100%): 


Pi 


In 


ftpHU,  50% 

1  {frpg-frp,  3) 

(ftp, 1~  ftp, 3) 


ft 


PHU,  50% 


+  P2-ftl 


PHU,  50% 


(74) 


P2  = 


In 


ft! 


1- 


PHU,  50% 

ftp 2  —  ftp, 3 


~ft! 


PHU,  50% 


ft P,l~  ft l 


P,3 


—  In 


1- 


ftpHU,  100% 
ftp,-  3  —  ftp, 3 

ftp,l-ftl 


-ft! 


PHU,  100% 


P,  3 


ftpHU,  100%  '  ftpHU,  50% 


(75) 


where  the  normal  fraction  of  phosphorus  in  the  plant  biomass  near 
maturity  is  used  in  Equation  75  to  ensure  that  the  denominator  term 

(  ftp,  ~3  ftp,  3 

(ftp,l-ftp,3] 

(frp,-3- ftp, s)  =0.00001. 


does  not  equal  1.  The  formulation  assumes 


To  determine  the  mass  of  phosphorus  that  should  be  stored  in  the  plant 
biomass  for  the  growth  stage,  the  phosphorus  fraction  is  multiplied  by  the 
total  plant  biomass: 


biop,oPt=ftpbio  (76) 

The  plant  nitrogen  demand  for  a  given  day  is  a  function  of  the  difference 
between  the  phosphorus  content  of  the  plant  biomass  expected  for  the 
plant’s  growth  stage  and  the  actual  phosphorus  content: 
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Pup  =  1.5 [bioP  opt  - bioP )  (77) 

where  Pup  is  the  potential  phosphorus  uptake  (kg  P/ha).  The  difference 
between  the  phosphorus  content  of  the  plant  biomass  expected  for  the 
plant’s  growth  stage  and  the  actual  phosphorus  content  is  multiplied  by  1.5 
to  simulate  luxury  phosphorus  uptake,  the  continued  uptake  beyond  what 
is  required  for  immediate  growth. 

The  depth  distribution  of  phosphorus  uptake  is  calculated  with  the 
equation: 


up,z 


Pup 

1-exp 

f  V 

6  z 

l-exp(-pp) 

P p  z 

* root  J 

(78) 


where: 

Pup,z  =  potential  phosphorus  uptake  from  the  soil  surface  to  depth  z 
(kg  P/ha) 

(),,  =  phosphorus  uptake  distribution  parameter. 

The  potential  phosphorus  uptake  for  a  soil  layer  is  calculated  by  solving 
Equation  78  for  the  depth  at  the  upper  and  lower  boundaries  of  the  soil 
layer  and  taking  the  difference: 

P  =  p  _  p  n 

up,Iy  up,zl  up,zu  ' 


where: 

P uP,iy  =  potential  phosphorus  uptake  for  layer  ly  (kg  P/ha) 

PuP,zi  =  potential  phosphorus  uptake  from  the  soil  surface  to  the  lower 
boundary  of  the  soil  layer  (kg  P/ha) 

PuP,zu  =  potential  phosphorus  uptake  from  the  soil  surface  to  the  upper 
boundary  of  the  soil  layer  (kg  P/ha). 

Root  density  is  greatest  near  the  surface,  and  phosphorus  uptake  in  the 
upper  portion  of  the  soil  will  be  greater  than  in  the  lower  portion.  The 
depth  distribution  of  phosphorus  uptake  is  controlled  by  (),,,  the  phos¬ 
phorus  uptake  distribution  parameter,  a  variable  users  are  allowed  to 
adjust. 
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Phosphorus  removed  from  the  soil  by  plants  is  taken  from  the  solution 
phosphorus  pool.  The  importance  of  the  phosphorus  uptake  distribution 
parameter  lies  in  its  control  over  the  maximum  amount  of  solution  P 
removed  from  the  upper  layers.  Because  the  top  to  mm  of  the  soil  profile 
interacts  with  surface  runoff,  the  phosphorus  uptake  distribution  param¬ 
eter  will  influence  the  amount  of  labile  phosphorus  available  for  transport 
in  surface  runoff.  This  formulation  allows  lower  layers  in  the  root  zone  to 
fully  compensate  for  lack  of  solution  P  in  the  upper  layers,  so  there  should 
not  be  significant  changes  in  phosphorus  stress  with  variation  in  the  value 
used  for  pp. 

The  actual  amount  of  phosphorus  removed  from  a  soil  layer  is  calculated: 

^ actualup,ly  min i^PUpiy  T  ^demand’ ^solution, Zy  )  (80 


where: 


P actuaiuP,iy  =  actual  phosphorus  uptake  for  layer  ly  (kg  P/ha) 

P demand  =  phosphorus  uptake  demand  not  met  by  overlying  soil  layers 
(kg  P/ha) 

P solution, iy  =  phosphorus  content  of  the  soil  solution  in  layer  ly  (kg  P/ha). 

Crop  yield  and  nutrients  removed 

When  a  harvest  or  harvest/kill  operation  is  performed,  a  portion  of  the 
plant  biomass  is  removed  from  the  land  surface  as  yield.  The  nutrients  and 
plant  material  contained  in  the  yield  are  assumed  to  be  lost  from  the 
system  (i.e.,  watershed)  and  will  not  be  added  to  residue  and  organic 
nutrient  pools  in  the  soil  with  the  remainder  of  the  plant  material.  In  con¬ 
trast,  a  kill  operation  converts  all  biomass  to  residue. 

Crop  yield 

The  crop  yield  is  calculated  as: 


yld  =  bioag  ■  HI 


yld  =  bio- 


1 +  HI) 


HI  <  1.0 


(81) 


HI  >  1.0 


(82) 
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where: 

yld  =  crop  yield  (kg/ha) 

bioag  =  aboveground  biomass  on  the  day  of  harvest  (kg/ha) 

HI  =  potential  harvest  index. 

The  aboveground  biomass  is  calculated: 

bi°ag  =  (1-frroot)bio  (S3) 

where  fr root  is  the  fraction  of  the  total  biomass  in  the  roots  the  day  of 
harvest. 

Potential  harvest  index 

The  fraction  of  the  aboveground  plant  dry  biomass  removed  as  dry 
economic  yield  is  called  the  harvest  index.  For  the  majority  of  crops,  the 
harvest  index  will  be  between  o.o  and  l.o.  However,  plants  whose  roots 
are  harvested,  such  as  sweet  potatoes,  may  have  a  harvest  index  greater 
than  1.0. 


The  economic  yield  of  most  commercial  crops  is  the  reproductive  portion 
of  the  plant.  Decades  of  crop  breeding  have  led  to  cultivars  and  hybrids 
having  maximized  harvest  indices.  Often,  the  harvest  index  is  relatively 
stable  across  a  range  of  environmental  conditions. 


The  harvest  index,  for  each  day  of  the  plant’s  growing  season,  uses  the 
following  equation: 


jjj  _  jjj  _ 100  frpHU _ 

100  frPHU  +  exp  (ll.l  — 10  frPHU) 


(84) 


where: 

HI  =  potential  harvest  index  for  a  given  day 
Hlopt  =  potential  harvest  index  for  the  plant  at  maturity  given  ideal 
growing  conditions. 
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Actual  harvest  index 

The  harvest  index  is  affected  by  water  deficit;  actual  harvest  index  is 
calculated  using  the  following  equation: 


HIact={HI-HImin) 


Yt 


Yu>u  +exP(6.13  — 0.833Yum) 


HT 


(85) 


where: 

Hlact  =  actual  harvest  index 

HImin  =  harvest  index  for  the  plant  in  drought  conditions  and 

represents  the  minimum  harvest  index  allowed  for  the  plant 
Yum  =  water  deficiency  factor. 

The  water  deficiency  factor  is  calculated: 

771 

E*. 

Y„=100*-b!—  (86) 

EE° 

i= 1 

where: 

m  =  day  of  harvest  if  the  plant  is  harvested  before  it  reaches 

maturity  or  the  last  day  of  the  growing  season  if  the  plant  is 
harvested  after  it  reaches  maturity 
Ea  =  actual  evapotranspiration 
E0  =  potential  evapotranspiration. 
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Nutrients  removed 

The  amount  of  nutrients  removed  in  the  yield  is  calculated  as: 


Nyld=frN,yld-yld 

(87) 

Pyld  =  fi'p,  yld  ■  lJld 

(88) 

where: 

Nyid  =  amount  of  nitrogen  removed  in  the  yield  (kg  N/ha) 
frN,yid  =  fraction  of  nitrogen  in  the  yield 

Pyid  =  amount  of  phosphorus  removed  in  the  yield  (kg  P/ha) 
frp,yid  =  fraction  of  phosphorus  in  the  yield  (kg/ha). 

If  the  harvest  index  override  is  used  in  the  harvest  only  operation,  the 
formulation  assumes  that  a  significant  portion  of  the  plant  biomass  is 
being  removed  in  addition  to  the  seed.  Therefore,  instead  of  using  the 
nitrogen  and  phosphorus  yield  fractions  from  the  plant  growth  database, 
the  formulation  uses  the  total  biomass  nitrogen  and  phosphorus  fractions 
to  determine  the  amount  of  nitrogen  and  phosphorus  removed: 


Nyid  =  K-yld 

(89) 

Pyid  =frP- yld 

(90) 
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6  Model  Validation  Studies 

In  developing  any  numerical  model,  validation  of  the  process  descriptions 
is  very  important  if  the  model  is  to  be  used  in  a  predictive  mode.  With  that 
in  mind,  two  demonstration  projects  are  currently  underway.  The  first  is 
the  Eau  Galle  watershed  located  in  Wisconsin  and  the  second  is  the  Upper 
Auglaize  Watershed  located  in  Ohio.  In  both  cases  model  results  were  not 
completed  in  time  for  this  publication.  This  report  will  give  a  general 
description  of  the  demonstration  sites  and  goals  of  the  studies.  Results  will 
be  published  later  in  companion  technical  notes. 

Eau  Galle  Watershed  demonstration 

The  Eau  Galle  River  watershed  encompasses  a  402-km2  area  in  northwest 
Wisconsin  (Figure  3  (Downer  2007)). 


Figure  3.  Location  of  Eau  Galle  River  watershed. 
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While  the  larger  basin  is  important,  the  lower  portion  of  the  basin, 

Figure  4,  is  relatively  data  poor.  The  data,  particularly  rainfall,  are  not 
sufficient  for  the  testing  of  new  model  developments.  The  upper  portion  of 
the  watershed,  that  portion  above  Spring  Valley  Dam,  Figure  5,  has  been 
the  subject  of  intensive  past  studies,  and  compared  to  the  lower  portion  of 
the  watershed,  is  relatively  data  rich.  This  103-km2  sub-watershed  was 
selected  to  demonstrate  the  NSM  model,  as  well  as  other  features  in  the 
GSSHA  model. 


Figure  4.  Spring  Valley  Dam  watershed. 
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In  addition,  the  Spring  Valley  Dam  watershed  has  these  important 
features: 

•  Reservoir  with  controlled  discharge 

•  Mixed  hydrology  basin  with  contributions  from  surface  and  subsurface 

•  Agricultural  basin  with  diffuse  nutrient  sources 

•  Dentritic  stream  network 

With  the  observed  data  and  complex  features  found  within  the  Upper  Eau 
Galle  Watershed,  this  study  area  is  ideal  for  testing,  validating,  and  further 
developing  the  flow,  sediment,  and  nutrient  algorithms  found  within 
GSSHA  and  NSM. 

Upper  Auglaize  Watershed  demonstration 

The  Upper  Auglaize  Watershed  is  located  in  portions  of  Auglaize,  Allen, 
Putnam,  and  VanWert  Counties  in  the  southern  portion  of  the  Maumee 
River  Basin.  The  watershed  encompasses  211,956  acres  (-330  mi2) 
upstream  of  the  USGS  Fort  Jennings  gauging  station,  Figure  5.  Land  use  is 
predominantly  agricultural  with  74.2  percent  cropland,  10.8  percent 
grassland,  6.2  percent  woodland,  and  8.8  percent  urban  and  other  land- 
uses.  Blount  and  Pewamo  are  major  soil  types  in  the  watershed.  These 
soils  are  characterized  as  somewhat  poorly  to  very  poorly  drained  with 
moderately  slow  permeability.  Farm  fields  in  the  watershed  are  extensively 
tile  drained  and  the  area  has  a  very  extensive  network  of  man-made  or 
man-altered  drainage  ditches.  Common  conservation  practices  applied  in 
the  watershed  include  grassed  waterways,  tile  and  surface  drainage, 
conservation-tillage  and  no-tillage,  grass  filter  strips,  and  erosion  control 
structures.  Soils  are  among  the  most  productive  in  Ohio  and  fields  not  still 
in  woodland  or  USDA  set-aside  cover  are  intensively  farmed.  At  the 
present  time,  no-tillage  (no-till)  is  practiced  on  51  percent  of  the  cropped 
fields  and  mulch  tillage  on  17  percent.  Most  natural  fencerows  have  been 
removed.  There  is  a  modest  amount  of  livestock  in  the  watershed,  most  of 
which  is  in  confined  feeding  setups  with  very  few  animals  on  pasture. 
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Upper  Auglaize  Watershed 
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Figure  5.  Upper  Auglaize  Watershed  (Bingner  et  al.  2005). 
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The  overall  purpose  of  the  Upper  Auglaize  Watershed  Demonstration  is  to 
perform  the  monitoring  and  modeling  design  needed  to  create  both:  l)  the 
linkage  between  the  established,  recognized  edge-of-field  tool 
DRAINMOD  and  established  subwatershed  and  watershed  models,  and 
2)  the  capability,  at  the  subwatershed  and  watershed  scale,  to  reliably 
predict  the  transport  of  agricultural  loadings  from  the  edge-of-field  to  both 
subwatershed  and  watershed  outlets.  This  missing  capability  may  contain 
the  understanding  of  system  response  required  to  effectively  predict,  track, 
and  manage  hypoxic  conditions  in  receiving  water  bodies. 

The  modeling  phase  of  the  project  will  involve  the  development  of 
enhanced,  integrated  modeling  components  for  surface/subsurface  trans¬ 
port  of  sediments  and  nutrients  from  agricultural  fields  to  the  outlet  of  the 
watershed.  This  component  of  the  study  is  envisioned  as  a  multi-year 
effort  (3  years)  and  features  the  integration  of  state-of-the-art  monitoring, 
modeling,  and  remote  sensing  components. 
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14.  ABSTRACT 


As  research  continues,  improved  understanding  of  nutrient  process  descriptions  and  new  capabilities  will 
be  developed  and  as  such  will  be  integrated  into  the  NSM.  Future  releases  of  NSM  will  include:  1)  improved 
nutrient  kinetics  for  channels,  2)  improved  soil  plant  dynamics  (EDYS  sub-module),  3)  inclusion  of  carbon  cycle 
modules  for  soils  and  overland  flow,  and  4)  improved  nitrogen  and  phosphorus  kinetics  for  soils  and  overland 
flow. 


